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RESEARCH SUMMARY 


Maximum load-loss rate within the combustion zone of a vertically 
(downward) spreading fire was obtained for excelsior (0.07 em in cross- 
section) at bulk densities from 0.0016 to 0.026 g/em®, Fuel was contained 
within a continuously weighed circular wire mesh basket 1 ft 2 (929 cm?) in 
circular area and 0.5 ft (15.2 cm) deep, surrounded by a separate annular 
basket of fuel, to minimize edge effects. Load-loss rates were similar to 
those obtained by Rothermel (1972) for horizontally spreading fires in the 
same fuel at the same bulk densities. 


The peak rate of load loss per unit, area exhibits a strong dependence 
upon bulk density for both vertically and horizontally spreading fires. The 
maximum of this peak rate occurs at the same bulk density in the two cases. 
Because of these similarities, this method offers an alternative for investi- 
gating the combustion dynamics of fire; moreover, the vertical array is 
much easier to use than the horizontal array. 


INTRODUCTION 


Frandsen and Rothermel (1972) have presented a method of measuring the energy- 
release rate of a spreading fire through the load-loss rate of a section of the fuel 
bed. This paper describes an alternative fuel bed that uses less fuel, that takes less 
time toyeonstruct, and that 1s less complex to measure than the fuel bed used previ- 
ously, Uel arrays to be examined are contained within a circular basket. The top is 
ignited uniformly over the surface and the basket's fuel load loss is continuously 
recorded as the fire moves vertically downward. 


Rothermel (1972) has shown that the fundamental dynamics of fire spreading through 
a porous fuel array can be expressed through T, the reaction velocity, the ratio of the 
efficiency of mass degradation during pyrolysis to the elapsed time of pyrolysis exper- 
ienced during fire spread. Its product with w_, the net original ovendry load of a fuel 
bed, is equivalent to w, the maximum load-loss rate of a portion of the fuel bed as the 
fire spreads through the fuel: 


w=wle 
n 


The method of obtaining W experimentally has been covered by Frandsen and Rothermel 
(1972) in their investigation of the reaction intensity of the combustion zone of a 
spreading fire. (The rate w is properly a loss; however, the negative sense is inter- 
nally compensated.) The reaction intensity, I,, is another important dynamic parameter 
Ofethe rire. It is the total energy-release rate evolved from the combustion zone, and 
is related to w through h, the low heat of combustion: 

In = hw = hw 
The load-loss rate, W; is an important basic parameter that relates to the dynamics of 
a spreading fire. The maximum load-loss rate of a burning fuel array can be obtained 
by monitoring the load of a basket of fuel as the fire moves with a horizontal plane 
downward through the basket (fig. 1). The purpose of this paper is to show that, for 
excelsior, the fuel basket experiences a maximum load-loss rate similar to that expe- 
rienced by a horizontal fuel bed slice Ax, (fig. 2) at the same packing ratio, the 
ratio of the bulk density to the particle density. Although the bed fire travels 
horizontally, the observer at Ax views a downward moving front that is similar to the 
front illustrated in figure 1; namely, as Ax tends to zero, the boundaries of the pyrol- 
ysis zone can be easily approximated as horizontal planes. Because only the fuel slice 
is weighed, the rate of change of the load within the slice is hypothesized to be 
identical to the rate of change of the load in the burning basket. 


Figure 1.--A horizontal pyrolysis 
zone moving downward at a rate, 
C, through a wetghed fuel basket 
of depth, &. A doughnut-shaped 
basket surrounds the wetghed 
basket to elimtnate edge effects. 


Figure 2.--A cross-sectional vtew of a fire spreading horizontally through a fine fuel 
(complete burnout). A wetghed fuel slice ts illustrated by Ax. The lower boundary 
of the pyrolysts zone ts the tgnitton tnterface with the fuel. 
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EXPERIMENTAL PROCEDURE 


The experiment was conducted in the combustion chamber of the Northern Forest Fire 
Laboratory. The combustion chamber is connected to an environmental conditioning facil- 
ity capable of maintaining humidity within +1.5 percent and temperature within +33° F 
(+0.5° C). The entire series of experimental burns was run at 80° F (27° C) and 20 
percent relative humidity, resulting in a 5 percent equilibrium moisture content for 
the fuel. The burns were conducted under an exhaust flue in the center of the combus- 
tion chamber 44 by 44 by 66 ft (13.4 by 13.4 by 20.1 m.) 


A wire mesh basket (13 by13 mm grid mesh size), 1 ft? (929 cm*) in circular area 
and 0.5 ft (15.2 cm) deep was used to contain the fuel. Square cut aspen wood excelsior 
(Populus tremulotdes Michx.), 0.07 cm in cross-sectional dimension (a fine fuel) was 
used; the same fuel was used in the original spreading fire experiment to which the 
burning basket results are to be compared. The particle density was 0.40 g/cm’. Dif- 
ferent packing ratios were obtained by changing the load while holding the depth fixed. 
A doughnut-shaped basket with an annular space of 0.5 ft (15.2 cm) and identical depth 
and packing ratio surrounded the inner basket (fig. 1) to eliminate edge effects and 
thus maintain a horizontal planar combustion front as the fire moves downward through 
the inner basket. Separation between the baskets is highly important. An annular space 
of one-fourth inch (0.6 cm) separating the inner and outer baskets was optimum for these 
tests. A larger space will result in increased burning within the open annular space 
deteriorating the horizontal combustion front, whereas a smaller space will result in 
serious disturbance to the load measurement through connecting bridges occasioned by 
heat warping the adjacent wire mesh walls of the baskets. 


A string presoaked in xylene, a highly flammable fluid, was laid back and fourth 
over the upper surface of the fuel to insure uniform ignition. The load was monitored 
by a load cell transducer located below the inner fuel basket and insulated from the 
heat by baffles that also prevented airflow through the annular space between the center 
basket and the outer ring of fuel. 


RESULTS 


x 


Figure 3 is a representative set of fuel load histories for each packing ratio, 8. 
The load-loss rate was obtained simultaneously with the load history by electronically 
differentiating the signal from the load cell transducer. The differentiated signal 
was also used to determine the first indication of weight loss for the curves in figure 
3 because of its greater sensitivity to changes in the load. Except for 8 = 0.065, each 
curve has an initial negligible load loss. 
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Figure 3.--Typtcal load histortes of burning fuel baskets at six different packing 
ratios. The percent load loss through combustion ts indtcated at the end of the 
first major pertod of fuel consumptton and at 30 seconds for 8 = 0.004 to 8 = 0.016 
and at 50 seconds for — = 0.032 to 8 = 0.068. 


Five replications of the load history were made at six packing ratios from 0.004 
to 0.065. Table 1 is a compilation of the maximum load-loss rate, w, obtained from 
the burning basket. Table 2 contains Rothermel's original data for excelsior from 
horizontal spreading fires. Both sets of data are displayed in figure 4 along with a 
prediction of w according to Rothermel (1972) (appendix B). The prediction differs 
from the original excelsior data because the results of other fuel sizes were amalga- 
mated by Rothermel to obtain the generalized expression for reaction velocity and thus 
the maximum load-loss rate. 
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Figure 4.--A plot of the maximum load-loss rate, w, versus the packing ratto, 8, for 
excelstor burned tn a basket and in a fuel bed. The continuous curve was obtained 
from the expression for w in appendix A. The upper and lower extreme are the standard 
devtatton from the mean. 


DISCUSSION 


The data contained in tables 1 and 2 for the two different methods of obtaining 
the maximum load-loss rate do not lend themselves to a simple comparison at the same 
packing ratios. However, assuming that a quadratic model fits these data, we can 
approximate each data set with a quadratic regression curve, w regressed on £8, and com- 
pare these curves to a quadratic regression of the composite of both data sets. The 
composite curve expresses fuel loss rate as accurately as the individual data sets if 
the composite curve has a sum of squared errors (SSE) (squared differences between 
observed and predicted from quadratic model) not significantly greater than the sum of 
SSE's from the separate curves. It then follows that the curves are not dissimilar. 

A statistical test of the data from tables 1 and 2 (appendix C) shows that it is very 
unlikely that the two data sets are different. 


The early load history (fig. 3) of the burning fuel follows a declining sigmoid 
curve. Shape of the curve varies with packing ratio and is distinguished from other 
curves through the maximum load-loss rate occurring at the inflection. It is important 
to determine if the depth was sufficient to achieve the characteristic maximum load-loss 
rate for each packing ratio. Fuel consumption efficiency is a means of examining this 


question. It is indicated in figure 3 at the end of the curve and in the region of 
increasing slope (second knee except at 8 = 0.65) following the inflection of the 
Imttiale sapmord curye. For 6 = 0.032 to B = 0.065, 39 to 79 percent of the fuel still 
remains at 50 seconds--well after the initial sigmoid load history. Consequently, the 
depth. 15.2 cm, was sufticrent for 6 = 0.032 to 6 = 0.065. For 6 = 0.004 to 8 = 0.016, 
as little as 4 percent of the fuel remains at 30 seconds; however, each curve is sigmoid 
and shows at least 17 percent more fuel consumed beyond the second knee, indicating that 
the initial sigmoid curve was not foreshortened by lack of fuel depth. 


In the spreading fire, the major rate of fuel consumption is near the front of the 
combustion zone--consistent with its early occurrence in the fuel basket--and is of 
prime importance to the propagation of the fire (Frandsen and Rothermel 1972). 


Converting an idealized load history to load-loss history provides us with an easy 
reference for comparing the dynamics of burning fuel baskets of constant 8 having 
different depths (initial loads). Because the fuel below the pyrolysis zone has no 
effect upon the present load-loss history, we can see that the load-loss curve for an 
extreme depth will have all the characteristics of load-loss curves of lesser depth 
(fig. 5), except that the period of constant load-loss rate becomes longer. The slope 
of the curves in figure 5 to the right of the second knee is dependent on char buildup. 
Char formation depends on the mineral content of the fuel and increases with oxygen 
depletion that in turn increases with time as the pyrolysis zone travels down into the 
basket. The char buildup for excelsior is negligible up to and including the maximum 
load-loss rate because the maximum rate occurs early in the load history and the mineral 
content 1s low (0.03 percent). Whatever the effect, it is the same for both the slice 
and the basket. 


Load loss 


Figure 5.--An idealized load- 
loss htstory of a burning 
fuel basket of extreme 
depth, curve D, having 
similar characteristics 
to a basket of lesser 
depth, curves C down to A. 


Time 


The statistical test in appendix C indicates that the alternative method presented 
here can be used for fuel near the size of excelsior to obtain the reaction velocity 
used by Rothermel to describe the dynamics of a horizontally spreading fire. © Investi- 


gations at larger fuel sizes are necessary to extend the range of the usefulness of 
this method. 


This method should be useful in studying fire spread in forest fuels comparable in 
size to excelsior. Some specific applications that come to mind are; dependence of fire 
spread on natural variables such as mineral content, moisture content, extractive 
content, and packing ratio, and also the efficacy of fire-retarding chemicals. 
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APPENDIX 


SYMBOL 


Appendix A 


Terms 


DEFINITION 
Maximum load-loss rate 
Net original ovendry load (less inorganics) 
Reaction velocity 
Optimum reaction velocity 
Maximum optimum reaction velocity 
Low heat content of fuel 
Reaction intensity 
Fuel depth 
Surface area-to-volume ratio 


Packing ratio, the fraction of bulk volume 
occupied by fuel 


Optimum packing ratio 
Ovendry particle density 
Mineral damping coefficient 
Moisture damping coefficient 
Silica free ash content 


Fuel moisture content 


10 


cal/g 
cal/cm2-s 


cm 


Fraction 
of ovendry 
weight 


Appendix B 
Derivation of Maximum Load-Loss Rate 


Beginning with w = wis the net load 1s expressed as 


Wn epee 
where 
Pp = ovendry particle density 
6 = depth of the fuel basket or fuel slice 
8 = packing ratio, the fraction of the bulk volume occupied by fuel. 


The other term in the expression for w is the reaction velocity, T. Rothermel (1972) 
provides a set of parametric equations for obtaining IT from initial fuel parameters. 


te non! 
where 
N= mineral damping coefficient = f (silica free ash content) 
ies moisture damping coefficient = f (fuel moisture content) 
[* = optimum reaction velocity 
and 
r- = T* (6/8 )Xexp[A(1-8/8 J] 
max op op 
Beets (405) 3-0,0504 og) - =) 2 
max : 
B Sods ome tee 
op : 
A =i. 774 o * ~ 7.27) 
0) = particle surface area-to-volume ratio 


Substituting Wa and T into w gives 
Wes pon m8 I 
From Rothermel's (1972) original data 


26.12 lbs/ft? 


D 
' 


p 

Gy =p OR S75 ft 

a ee 0.694 @ Me = 0.05 

Nn. = 0.507 @S. = 0.0036 
Ss e 


Substituting into w we have 

w = 0.0280 8 I’ 
with o = 1,848 ft ! for excelsior 

P* = 15.237 e0- 3498 [g- e 8 10.3498 
where 


B~ = B/0.00707 
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Appendix C 


Statistical Test 


Statistical test of the similarity between the maximum 
load-loss rate dependence on packing ratio of the alter- 
native burning basket method and the horizontally 
spreading fire. 


The SSE and degrees of freedom, d.f., are given in the 
following tabulation. The model assumes a second degree 


polynomial. 
Data set SSE* dit. y=" N=s, 
Verticle burn, SSE(V) 9.5600x10 © 27, 
Horizontal burn, SSE(H) 43.6092x10 © 20 
Composite, SSE(C) 58.0511x10 & 50 


*Sum of squared errors 
Forming F value according to Draper and Smith (1966) 


(ASSE/Ad.f.) 


F © (SSE(V+H) /d.£. (V+H)) 
Therefore 
ar (4.8819x10 ©/3) _ (4.8819) (47) 
(53.1692x10 ©/47) (53.1692) (3) 
F = 1.4385 


Given no difference in the two sets of data the probability of such a value for F, 
is 0.24. Therefore, we accept the hypothesis that both methods are estimating the same 
result. 
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